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^> ' Abstract. Using Monte Carlo methods to simulate the inverse Compton scatter- 

. ing of soft photons, we model the spectrum of the Galactic black hole candidate 

Cyg X-1, which shows evidence of a nonthermal tail extending beyond a few 
. hundred keV. We assume an ad hoc sphere of leptons, whose energy distribu- 

' tion consists of a Maxwellian plus a high energy power-law tail, and inject 0.5 

. keV blackbody photons. The spectral data is used to constrain the nonthermal 

plasma fraction and the power- law index assuming a reasonable Maxwellian tem- 
perature and Thomson depth. A small but non-negligible fraction of nonthermal 
Oh' leptons is needed to explain the power-law tail. 

o 

^ : INTRODUCTION 



X 



Cygnus X-1 is the brightest and most studied Galactic black hole candidate. 
Its high-energy spectrum has been extensively modeled in attempts to deter- 
\ mine what processes produce the continuum emission. While several models 

debate the detailed geometry of the system [1-4], many agree that the 10- 
200 keV portion of the spectrum is a result of inverse Compton scattering of 
soft X-ray photons. Early modeling attempts assumed a single temperature 
Maxwellian plasma [5] . However such simple models do not adequately fit the 
observations [6,7]. More complex models, such as one with two plasma regions 
of different temperature, increase the goodness-of-fit substantially [3]. 

When Cyg X-1 is in the normal (72) state, a high-energy power-law is ob- 
served in the 7-ray continuum [8]. This is typically attributed to some non- 
thermal process. It is possible that the power-law is due to inverse Compton 
scattering off of a high-energy power-law tail in the plasma energy distribu- 
tion. If so, both the emission above and below 300 keV could originate from 
the same plasma. 



FITTING THE CYG X-1 SPECTRUM 



We wish to see if the high-energy spectrum of Cyg X-1 can be explained sim- 
ply by Comptonization from a single plasma region with a Maxwellian-hpower- 
law energy distribution. We simulated Comptonized spectra with a Monte 
Carlo code based on the algorithm of Pozdnyakov, Sobol, and Syunyaev [9]. 
We compared the results to combined BATSE and COMPTEL data from 
1991 [8,10]. Though this data has already been unfolded through the detector 
response using an assumed model, it serves to give us an approximate solution. 

The parameters that define the shapes of our simulated spectra are the 
temperature of the thermal leptons (electrons and pairs) Te, the Thomson 
depth Tt, the fraction of nonthermal leptons ^, and the energy index of the 
nonthermal leptons p. Generating and interpolating over a 4-dimensional pa- 
rameter space is computationally intensive and for this study we fixed both 
Te = 65 keV and tt = 2.45. A value of Te very close to this was deter- 
mined for this data assuming single temperature analytic Comptonization [8]. 
The parameter tt was determined from the spectral index a. of the 30 to 70 
keV spectra, as suggested by Pozdnyakov, Sobol and Syunyaev [9], where the 
equations 



are solved for tt. This left two free spectral shape parameters, ^ and p. We 
calculated a 6 x 6 grid of simulated spectra, with ^ ranging logarithmically 
from 0.25% to 8.0% and p ranging from 3.25 to 4.5. These grid points were 
chosen based on our experience with this code applied to other astrophysical 
phenomena [11,12]. We iterated our code until the statistical signal-to-noise 
within the range of 30 keV to 2020 keV was less than 10 in each of the 22 bins. 
Bins above 500 keV were smoothed in a manner similar to that in Pozdnyakov, 
Sobol and Syunyaev [9]. A simple spline algorithm allowed us make our model 
continuous so we could evaluate a between our simulation and the discrete 
data. 

Figure 1 shows a reasonable fit of our code to the BATSE and COMPTEL 
data, where ^ = 0.5% and p = 3.5. To see what range of the parameter space 
is acceptable, we next examined over the entire grid. In Figure 2, we plot 
the confidence levels of fits on this grid. Contours are drawn here assuming 
a polynomial interpolation between grid points.. This shows that for a fixed 
Tc = 65 keV and tt = 2.45, there is a 68.3% (la) confidence that ^ lies 
between 0.25% and 1.0% and p is between 3.25 and 4.25. 
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FIGURE 1. Fit of Maxwellian+nonthermal Comptonization model to unfolded spectrum 
of Cygnus X-1 where To = 65 keV, tt = 2.45, C = 0.5%, and p 3.5. Also plotted for 
comparison is a fit of the same model with no nonthermal electrons or pairs. 



We remind the reader that this is a crude way to test a model since we are 
comparing to aheady unfolded data. However, by showing a reasonable match 
between our simulated spectra and the unfolded spectra, we can be confident 
that more a rigorous procedure would also work. 



SUMMARY AND CAVEATS 

We find that the X-ray /7-ray spectra of the normal state of Cygnus X-1 
can roughly be reproduced by Comptonization of 0.5 keV l)lackbody photons 
through a combination thermal-nonthermal plasma. For Tc=65 keV and the 
Thomson depth tt=2A5, we determine that the fraction of nonthermal leptons 
^ = 0.5%^o;25% ^^"^ the energy index of the nonthermal leptons p = S.5li.oll- 

While these results place first-order hmits on the nonthermal lepton distri- 
bution, several modifications to these procedures would be necessary in order 
to reliably determine the allowed parameter space. The current Monte Carlo 
code does not require self-consistency between the lepton distribution and the 
radiative cooling. This obviously must be corrected to produce physically 
meaningful results. It is also necessary to generate a four-dimensional grid 
of simulated spectra to allow all four shape parameters to vary. This grid 
should also extend to lower values of ^ and p since our current grid does nec- 
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FIGURE 2. Confidence level contours in ^ — p parameter space for fits of model to Cyg 
X-1 spectrum. Confidence levels are calculated from goodness-of-fit statistic assuming 
11 degrees of freedom (14 data points - 3 fit parameters). See Press et al. 1993 for details. 

essarily exclude solutions in this portion of the parameter space. Finally, any 
future work would require folding our simulated spectra through the detector 
responses to allow direct comparison with the 7-ray count data. 
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